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Oblique incidence of a p-polarized laser beam on a fully ionized plasma with a low density plasma 
corona is investigated numerically by Particle-In-Cell and Vlasov simulations in two dimensions. A 
single narrow self-focused current jet of energetic electrons is observed to be projected into the 
corona nearly normal to the target. Magnetic fields enhance the penetration depth of the electrons 
into the corona. A scaling law for the angle of the ejected electrons with incident laser intensity is 
given. 
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The availability of tabletop high intensity laser systems has lead to the investigation of novel regimes of short pulse 
laser-plasma interaction. Recently the emission of collimated electron jets under specular angles with respect to the 
density normal direction have been observed for an obliquely incident laser beam on a steep density plasma 0. 

When a target is irradiated by an intense laser pulse above the field ionization threshold it quickly ionizes [0. 
For sufficiently long laser pulse irradiation the plasma present on the surface gradually expands into the vacuum 
with the ion acoustic speed. Hence, a plasma corona is formed. For short laser pulses there is not enough time for 
hydrodynamic expansion. Short pulse simulations show however, that an ion shelf is formed on a typical time scale 
t s = Lu ~ 1 (mi/Zm c ) 1 / 2 due to the generation of a strong electric field at the plasma-vacuum interface [||. This ion 
shelf represents a low density plasma corona. 

There are different mechanisms which can lead to collimated electron jets when an intense laser pulse interacts 
with a vastly over-dense, steep density plasma that has a low density plasma corona. One effect that plays a role 
in the interaction is the Brunei effect |I| which works for oblique incidence. Here the electrons are accelerated into 
the vacuum as well as into the target by the electric field present along the density gradient. The coronal plasma is 
expected to collimate and enhance the range of the ejected electrons in two (2D) or three (3D) spatial dimensions by 
quasi-steady magnetic field generation. The collimating effect of quasi-steady fields has recently been addressed in 
HU in a different context. 

To investigate the phenomenon just outlined in more detail we perform Particle-In-Cell (PIC) and Vlasov (VL) 
simulations both in two spatial dimensions (2D). In both approaches we do not simulate the evolution of the corona 
self-consistently but treat it parametrically with ions fixed. In our PIC simulations we rotate the target and in our 
VL simulations we boost to the frame of normal incidence to model oblique incidence of the laser beam. The boost 
frame method is well established in ID [Q-p). However, in 2D using the boost frame is very helpful in establishing the 
physics of the underlying laser-plasma interaction. 

We investigate the interaction of a p-polarized laser beam incident under angles of 30° (PIC) and 45° (VL) on a 
preformed fully ionized target with an under-dense plasma corona in front of it. In both simulations the laser beam 
has a duration of about lOOfs. For the PIC case the laser beam intensity is 2.0 • 10 18 W/cm 2 and for the VL case 
1.0T0 17 W/cm 2 . The laser wavelength in the simulations is l^m with beam diameters of 8 /im (PIC) and 5 /xm (VL) at 
full-width-half-maximum. The coordinates of the simulation box are x and y, respectively. The size of the simulation 
box is 23 /im x 23 /im for the PIC simulations and 6 /im x 12 /im for the VL simulations. In our PIC simulations 
we assume mirror- reflecting boundary conditions for the electrons in x- and y-directions. In the VL simulations we 
use periodic boundary conditions in y-direction. Electrons leaving the simulation box at x = 6fim are replaced by 
thermal electrons and at x = 0/iin they are allowed to escape. We note that mirror-reflecting boundary conditions 
in our PIC simulations force us to increase the simulation box for long simulation times. The distribution functions 
for the electrons and ions needed for the VL simulations have two momentum directions p x and p y in addition. The 
quasi-particle number per cell used in the PIC simulations is 50 for each species. The fully ionized plasma density is 
4n c (PIC) and 8n c (VL). In both simulations we assume a uniform low density plasma corona with a density of 0.1n c 
in front of the target. 

Plot (a) of Figure [l] gives the quasi-steady magnetic field B z in front of the target inclined at 30° , as obtained by 
PIC simulations. The peak magnitude of the normalized magnetic field is 0.62, which corresponds to approximately 
30 MG. It changes polarity very rapidly along the density gradient revealing the presence of a very localized self- 
focused current jet. The low density plasma corona guarantees quasi-neutrality and helps to generate the magnetic 
field in front of the target. Plot (b) of the same figure shows the electron energy density. We find a collimated electron 
jet which coincides with the quasi-steady magnetic field from plot (a). For the parameters of plot (b) the ejection 
angle is approximately 17° form the target normal. There are also fast electrons injected into the over-dense plasma. 
We again observe that they are almost normal to the target surface. Figure |l| (c) shows the instantaneous plot of the 
electron energy density with over-plotted positive B z field indicating the phase of the laser field. It is clearly seen that 
the outgoing electrons are generated on the target surface once per laser cycle by the Brunei absorption mechanism 
||,|| and are bunched on the scale of the laser wavelength consequently. The range of the electrons is enhanced. A 
similar result we obtain from our VL simulations which make use of boost frame coordinates. 

To illustrate how the boost frame approach for oblique incidence in 2D works we briefly derive the correct boundary 
conditions for the laser pulse in the boosted frame. We start by defining an arbitrary pulse envelope function z(x, y, t) 
in the lab-frame. Next we perform a Lorentz rotation of electromagnetic fields about (xa,ya). In the final step we 
boost the latter to the frame of normal incidence for which the longitudinal field E x disappears. We obtain 

E*=0, Ey = i z (x r ,y r ,t) , B? = -4 z (x r ,y r ,t) , (1) 

where 



2 



~(x-x ) + (y-yo)P , y r = =(y-yo)-(x-x o )0, (2) 
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with 

i = ^V, x = -ct B , y=- 1V B . ,:-!) 
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The function z is the same function as in the lab- frame. For the relativistic factors we have (5 = sin and 7=1/ cos 0, 
where is the angle of incidence. Plot (a) of Figure |^ illustrates the incident time resolved electromagnetic field E y 
for a Gaussian pulse envelope. Plot (b) of the same figure gives the incident time resolved electromagnetic field E y of 
the simulations. 

Plot (a) of Figure |^ gives the quasi-steady magnetic field in the plasma corona in front of the over-dense plasma 
target. Plot (b) of the same figure gives the quasi-steady magnetic field with the quasi-steady B% over-plotted (red 
solid lines). Plot (c) of Figure || gives the quasi-steady magnetic field with the quasi-steady longitudinal current 
density j xe over-plotted (red dashed lines). 

Since the current density j xc is invariant under Lorentz boosts along y it may serve as a quantity from which to 
determine the direction of the electron jets. We now introduce the coordinates x = xB an d £ = y B + /3ct B which move 
along with the background plasma current present in the boosted frame. Since the time-averaged current density 
(j'xe) m the co-moving coordinates varies slowly with time we obtain (j B e ) (x B , y B , t B ) = (j B e ) (x>0- This yields 

<j£>(X,O = <J»e>(X,70 ■ (4) 

The direction of the collimated electron jets in the lab frame can now be calculated from the direction of the current 
density in the boosted frame. Plot (c) of Figure || gives (j B e )- The direction of the emitted electrons is indicated by 
the white solid line plotted in the figure. We obtain a mean emission angle of 20° in the boosted frame and 14° in 
the lab frame. We note that the lab frame is dilated in transverse direction when viewed from the boosted frame and 
hence the emission angle in the boost frame is larger by a factor of 7 = l/\/l — /3 2 as indicated by Equation (^). 

In boost frame coordinates we may easily analyze the physical mechanism that leads to the large areal quasi-steady 
magnetic field and the direction of the ejected electrons. We recall that in the boosted frame we have a constant 
background fluid velocity ub — csm9 which approaches speed of light for large angles of incidence. In this frame the 
polarization of the magnetic field vector of the incident laser beam is normal to the xy-plane and to the flow direction 
of the background current. If the laser intensity is small enough as in |jj and the angle of incidence sufficiently large 
the boost velocity exceeds the laser quiver velocity. The driving force under these conditions is exerted predominantly 
by the oscillating magnetic field of the laser beam (see the red solid contour lines of £> 2 plotted over B z in plot (b) 
of Figure |^ for the location of the force). The resulting force is F = — e ub x B and is capable of ejecting electrons 
out of the surface at a rate of once per laser cycle. This is the Brunei mechanism |Q. The quasi-steady magnetic field 
in the plasma corona is generated by the electron current emitted from the target. The polarization of the magnetic 
field is such that it collimates the electrons propagating through the plasma corona. 

To derive an approximate criterion for the angle range under which the fast electrons are emitted from the target 
surface we assume that the laser target interaction in the boosted frame is quasi-one-dimensional. Since the full- 
width-half-maximum of the laser beams in our simulations is at least 5 /im and the intensities are sufficiently low to 
prevent target imprinting we believe that this assumption is justified. We next rewrite the Vlasov equation in the 
boosted frame and solve it for an initial Maxwellian. We approximate the plasma-vacuum interface by a step-like 
density profile with nix) = uq for x > and treat the ions as immobile. We obtain for the distribution function 

__( P 2 x {0)+P 2 M\ ( ( Py (0) + pjmcf 



and for the equations of motion 



f(t) = p exp I exp - ^"V , (5) 



x(t) = x - / d-q v x {rj) , (6) 



Px(T)=p x + ej drj [E x ix(r]),ri) +v v ir)) d x A v ix(r]),ri)] , (7) 

P y ( T ) =P y + e [Ay(x(T), t) - A y (x, t)] , (8) 
p z {r)= Pz , (9) 
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with 



Vx/y(T) = . = • (10) 

ym 2 c 2 +pI(t) +p 2 y {T) +P 2 z (t) 

Equations (@) and (^) indicate lateral canonical momentum conservation in boost frame coordinates. We now 
assume that A y has a harmonic time dependence. Making use of Equations (||) and (||) and assuming « c or 
n x « c we obtain (p y ) s» —fijmc. The quantity (p y ) denotes the ensemble and time averaged transverse momentum. 
Treating (p x ) as a free parameter and transforming back to the lab frame yields 




The ejection angle is now given by tan# = (Py)/(p^)- For (p x ) -^oowe obtain tan# = 0j = tan#. This means 
that only ultrarelativistic electrons are ejected at very close to specular direction. For smaller longitudinal momenta 
(Px) we expect that the electrons are emitted at angles that are smaller than the angle for specular emission as observed 
in our simulations. Assuming that the mean fast electron momentum in rc-direction is given by (p x ) /jmc » V a IX 2 
we thus obtain 



, VI + al\ 2 - 1 

tanfl = ■ tan 6. (12) 

y/^TX 2 

Equation ( [l2|) looses validity as soon as target deformations start to become significant. The validity also depends 
on the accuracy of the mean longitudinal momentum given as a function of intensity. For IX 2 — 1.0 • 10 17 Wcm _2 /mi 2 
we obtain an ejection angle of 9 = 14° and for IX 2 = 2.0 • 10 18 Wcm~ 2 /im 2 wc obtain 8 = 17° from the simulations. 
This yields a" 1 8.0 • 10 17 Wcm- 2 ^m 2 . 

In conclusion, we have demonstrated with the help of two different simulation techniques that collimated electrons 
with enhanced range can be emitted from an over-dense target if a low density plasma corona is present. In addition, 
we have shown that fast electrons are injected into the over-dense plasma. Both, the ejection and injection directions 
are almost along the density normal direction for p-polarized light. By a transformation to the moving frame in which 
the laser pulse appears to be normally incident we were able to give a criterion for the angle range of the emitted 
electrons with ejection momentum. We find that for a planar interaction interface only speed of light electrons can 
be emitted at specular direction for p-polarized light. Less energetic electrons appear under almost normal emission 
angles due to a lack of lateral momentum transfer. This analytical result is in qualitative agreement with our numerical 
observations. We note that in addition to the mechanism outlined in this paper other mechanisms of fast electron 
generation like wake-field acceleration in the corona may exist leading to different emission angles. 
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FIG. 1. Quasi-steady magnetic field B z (a), quasi-steady electron energy density e (b), and instantaneous electron energy 
density with the laser field B z (c). Yellow contour areas are positive and blue areas negative. The parameters are n/n cr i t = 4, 
IX 2 = 2.0 • 10 18 Wcm _ Vm 2 , 6 = 30°, t = 120ft, and B z0 = 100MG. 
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FIG. 2. Illustration of the boost technique in 2D for an incident laser pulse of Gaussian pulse envelope (a) and time 
resolved E y taken in the simulations (b). The left figure in (a) gives the incident pulse, the figure in the middle gives the 
pulse after a rotation, and the right figure shows the pulse after the final boost. The arrows in (a) indicate the propagation 
direction of the laser pulse. The white solid lines in (b) give the density profile. The parameters for (b) are n/n cr i t = 8, 
IX 2 = 1.0 ■ I0 17 Wcm-Vni 2 , = 45°, and t = 25fs. 

FIG. 3. Quasi-steady B z (a), quasi-steady B z with quasi-steady B 2 over-plotted (b), and quasi-steady B z with current 
density over-plotted (c). Yellow contour areas are positive and blue areas negative. The white solid lines along y located 
between x = 4/nm and x = 5/nm are density iso-contours. The other lines over-plotted in (a) indicate the quasi-steady magnetic 
field. They are along x = 3.63/im (solid) and y — 7.18/nm (dashed). The parameters are n/n cr i t = 8, IX 2 = 1.0-10 17 Wcm -2 /im 2 , 
9 = 45°, t = 75fs, and B = I.5MG. 
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